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ABSTRACT 

Observations of surface magnetic fields are now within reach for many stellar types thanks to the 
development of Zeeman-Doppler Imaging. These observations are extremely useful for constraining 
rotational evolution models of stars, as well as for characterizing the generation of magnetic field. We 
recently demonstrated that the impact of coronal magnetic field topology on the rotational braking 
of a star can be parametrized with a scalar parameter: the open magnetic flux. However, without 
running costly numerical simulations of the stellar wind, reconstructing the coronal structure of the 
large scale magnetic field is not trivial. An alternative -broadly used in solar physics- is to extrapolate 
the surface magnetic field assuming a potential field in the corona, to describe the opening of the field 
lines by the magnetized wind. This technique relies on the definition of a so-called source surface 
radius, which is often fixed to the canonical value of 2.5Rq. However this value likely varies from star 
to star. To resolve this issue, we use our extended set of 2.5D wind simulations published in 2015, 
to provide a criteria for the opening of field lines as well as a simple tool to assess the source surface 
radius and the open magnetic flux. This allows us to derive the magnetic torque applied to the star 
by the wind from any spectropolarimetric observation. We conclude by discussing some estimations 
of spin-down time scales made using our technique, and compare them to observational requirements. 


1. INTRODUCTION 

The magnetic fields of most stars are created by con¬ 
vective motions and large scale flows in their envelopes, 
the source of a dynamo effect. Stellar parameters such 
as rotational period, mass and age influence this mag¬ 
netic activity. These internal processes are difficult to 
probe, and must be investigated through indirect tech¬ 
niques such as asteroseismology. However the surface 
manifestation of this magnetic field largely shapes the 
structure of stellar coronae. Studying the Sun’s corona 
greatly improved understanding of stellar atmospheres. 
Magnetic processes are thought to heat the corona up 
to several million Kelvin, hence driving a mag netized 
outflow into interplanetary space (Parker 1958). As a 
consequence, the Sun has an expanding atmosphere, the 
solar wind. Within the corona, the competition between 
the expanding outflow and magnetic forces leads to open 
field regions and closed magnetic loops or streamers. To 
reproduce this structure, models such as the potential 
field source surface model (PFSS) ( [Schatten et al.|[T9^ | 
have been developed. This model assumes a curreht tree 
magnetic field up to a source surface beyond which all 
the field lines are opened by the wind. Solar wind prop¬ 
erties observed by spacecraft at I AU, have lead to the 
development of empirical models using t he PFSS, such 
as the W ang-Sheeley-Arge (WSA) model (Wang & Shee- 
ley 1995), while efforts towards more self-consistent so¬ 
lar wind models have be en undertaken (see the review of 
Hansteen & Velli||2012|). The WSA model uses a value 
of 2.5 Kq as the radius of a spherical source surface 


(tss), which has been chosen to best match the polar- 


itv of the interplanetar's 

magnetic field observed at 1 AU 

(Altschuler & Newkirkj 

1969 

Hoeksema et al.||1983|) and 

has been extensiveR used ever since (Wang &: Sheele-y 

1994 

Schrijver & De Rosa 

2003 DeRosa et al. 

2012). 


However it has been proposed that this value should 


change during the Solar cycle (Lee et al. 2011) due to 
variations of the solar magnetic activity. 

Our knowledge of the Sun gives us precious insights 
into the coronae of other stars. Zeeman-Dopp ler Imag¬ 
ing ( Donati fc Brown||1997 Donati et al.||2006 ) uses the 
polarization of light in a line prohle, produced by Zee- 
man splitting, to study a stellar magnetic field. The ro- 
tationally modulated variability of that line profile, pro¬ 
vides information about the strength and geometry of 
the magnetic field. The deduced surface magnetic field 
can in turn be used as an input for coronal models, from 
which integrated parameters can be estimated. In partic¬ 
ular, observationally calibrated mass loss rates and open 
magnetic fluxes are key for better underst anding the ro¬ 
tation al evolution of stars. For instance, [Reville et al.| 
(2015) showed that the open flux is the relevant param- 


ing induced b'v stellar winds (ISchatzman 

19621 

Weber & 

Davis|1967 

Kawalei'l 

1988 

Matt et ai.|20l 

2). A 

full mag- 

netohydroc 

ynamic ( 

MHD) simulation is 

able to recover 


the coronal structure of such stars, however the much 
simpler PFSS model is likely to reproduce most of the 


large scale coronal magnetic field properties (Riley et al. 


2006). For this latter technique, the relevant parameter 
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-in addition to the surface field- is the source surface ra¬ 
dius. This model has been applied to ZDI targets, with 
source surface radii set to different values, sometimes 
thanks to promine nces observations, sometim es in a more 
arbitrary fashion (|Jardine et al.||2002| |2013|) . Given the 
differences of the stellar parameters: coronal tempera¬ 
ture, rotation rate, magnetic field strength and topology, 
is the fiducial solar value of 2.5 Rq (or 2.5i?* in a stel¬ 
lar context) a good choice ? How can one a priori set a 
reasonable value for the source surface radius ? 

We investigate this by comparing the 6 0 MHD sim¬ 
ulations performed in Reville et al.| (2015) to potential 
extrapolations and estimate the optimal source surface 
radii matching the open flux of the stellar coronae. This 
is described in Section In Section we propose a gen¬ 
eral method to estimate a priori the optimal rgs from 
stellar parameters, without running any simulations. We 
find that for fast rotators, magneto-centrifugal accelera¬ 
tion is key for assessing a corr ect value. We use a proce- 


( 1985| ) (detailed in appendix 0 


dure based on Sakurai 
to obtain the right velocity profile, taking into account 
rotation and magnetic field. In Sectionj^we discuss some 
applications of our method on spin-down time scale for 
young stars, and summarize our conclusions. 

2. COMPARISON OF SELF-CONSISTENT MHD 
WIND SIMULATIONS AND POTENTIAL 
MAGNETIC FIELD EXTRAPOLATION 

2.1. Wind simulations with the MHD code PLUTO 

For more than two decades MHD simulations have 
been used to study the properties of stellar winds 


( Washimi fc Shibata||I993| iKeppens fc Goe dbloed||1999 
Matt fe Halick||2UIJ4 fAIatt fe Pudritz||2L)08p . Computing 
■ has allowed tor the inclusion of complex magnetic 


power 


mensions (Cohen et al. 2011 Vidotto et al. 

2014 

Stru- 

garek et al.|j2U14). In 

Reville et al. 

(2U15|), we presented 


pact of the magnetic field topology on the stellar wind 
brakin g. We used the set of parameters of |Matt et al.| 
(2012) but extended it to more complex topologies than 
the dipole, such as the quadrupole and the octupole, as 
well as combinations of multipoles. From this study, we 
generalized the law giving the magnetic torque created 
by a wind on a solar-like star: 


ofIReville et al. I 

2015 

), computed with the PLUTO code 

(iVIignone et al. 

2007 

). All details about the numerical 


aspects of the study are given in Reville et al. (2015) 


especially the necessary boundary conditions to properly 
compute the torque created by the wind. However, the 
formulation 0 becomes useful to compute the torque of 
a given star only if $open is known from stellar param¬ 
eters. Running MHD simulations gives this value and a 
measure of the angular momentum loss. However, the 
goals here is to provide a simple method to compute this 
quantity without having to run time consuming simula¬ 
tions. 

A general method that has been used widely in the 
solar physics community is the potential field extrapola¬ 
tion, which recovers the structure of the magnetic field 
up to a source surface radius, and assumes the wind has 
made the field completely radial beyond this point. 

2.2. Potential extrapolation with a source surface 


Introduced by Schatten et al. (1969), the Potential 
Field Source Surface (Ph'SS) modk is able to extrapo¬ 
late the whole spatial structure of a magnetic held in a 
corona given the surface magnetic held. This model as¬ 
sumes that the magnetic held is current free in a shellular 
volume delimited by the stellar surface (of radius r = r*) 
and a source surface of radius Beyond this surface, 
the model mimics the effect of the wind, which opens held 
lines, by setting the magnetic held to be purely radial. 
Thus, in the region r* < r < rgg, we have: 


V X B = 0, 


( 3 ) 


hence there exists a scalar held $, a potential of the 
magnetic held which satishes: 


- V$ = B. 


( 4 ) 


Since V • B = 0, is a solution of the Laplace’s equa¬ 
tion and we can write: 


A$ = 0. 


( 5 ) 


Two conditions hx the value of the potential and thus 
the magnetic held in the whole domain. First the poten¬ 
tial should match the observed held at the surface of the 
star: 


= M. 


l-2m 




,2-4m^2 


Ueg.(l + /V^|)l/2 I 


^2 

open 


( 1 ) 


where Mw is the mass loss due to the wind, i?*, D*, and 
/ = are the stellar radius, rotation 

rate and break-up ratio, respectively, it's, K 4 and m 
are the htted parameters for the braking law. ^open is 
the value of the unsigned magnetic hux if the integra¬ 
tion surface Sr contains all closed magnetic loops (and is 
therefore a constant): 


<i>(r) = [ \B- dS\ (2) 

JSr 

This formulation and the associated coefhcients have 
been derived using the grid of 60 numerical simulations 


-^\r=r, = -Br{r^,9,(l)). ( 6 ) 

Then the potential must not depend on (0, ^) at r = 
Tss since the held is purely radial from this point. This 
implies 


^’(?' = Tss) = 0 . 


( 7 ) 


We used the derivation given in |Schrijver fc De Rosa| 
(2003) to implement our reconstruction. In this model 
the niagnetic held only depends on the stellar surface 
held and on the value of the source surface radius r^g. 
Beyond this point the wind has been able to open all the 
held lines, which assumes that the thermal and turbulent 
pressure are enough to counter the tension of the coronal 
magnetic held. We will see in the next section how to 
predict a value for rgg. 
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Fig. 1.— Comparison of the magnetic field lines obtained by the simulations in white and by the potential extrapolation with the optimal 
source surface (cyan lines), in the dipolar case (left panel) and the quadrupolar case (right panel). The magnetic loops are well reproduced, 
but the field lines obtained by the simulations are not purely radial right beyond rss,opt whereas they are -by construction- with the 
potential extrapolation. The cyan dashed line represent the optimal spherical source surface. Color background is the logarithm of the 
density. Grey and white lines are the sonic and the Alfven surfaces respectively, note that the Alfven surface is always beyond the source 
surface. 




2.3. Criteria for an optimal source surfaee 

In order to assess what is the best value for and try 
to predict it, we need to have a reference value. For this 
we define an “optimal” source surface that best describes 
the results of our 60 simulations. Since the quantity we 
are interested in is the open flux, we call the optimal 
source surface radius the zero of the function: 


^{rss') — ^ open{r ss') ^open,sim ( 8 ) 

There is a unique solution since the open flux ob¬ 
tained through the potential extrapolation is a decreas¬ 
ing function of that starts at the surface flux value 
at i^ss = and tends to zero as Cgg tends to infinity. 
We have been through all our simulations, and Table 

gives the corresponding optimal source surface radius 
as a function of the stellar parameters. The parameter 
VAlvesc = B^j^/4TFpfj^/2GM^Jr\ is the Alfven speed 
on the surface at the equator over the escape velocity, 
which characterizes the magnetic field strength, while 
the break-up ratio / characterizes the rotation rate. The 
other parameters used for the simulations have been kept 
fixed, the sound speed at the base of the corona over 
the escape velocity has the value Cg/uggc = 0.222 and 
7 = 1.05. The optimal rgg can be easily found thanks to 
a bisection or a numerical Newton-Raphson method. 

Figure U shows a comparison of the magnetic field ob¬ 
tained with the simulation (white lines) and with the po¬ 
tential extrapolation (cyan lines) using the optimal rgg 
for dipolar and quadrupolar topologies. It can be seen 
that the closed magnetic loops are well reproduced. The 
location of the optimal source surface fairly well matches 
the size of the largest closed coronal loop of the simu¬ 
lation even though this is not the criteria we chose to 
define it. 

A difference can be noticed when looking at open field 
lines, while the potential extrapolation of the magnetic 


TABLE 1 

Table of parameters and computed optimal 

SOURCE SURFACE 


Case 

'VaI'^csc 

/ 


^ ss,opt 





Dip. 

Quad. 

Oct. 

1 

0.0753 

9.95e-5 

5.2 

3.3 

2.7 

2 

0.301 

9.95e-5 

8.6 

4.7 

3.5 

3 

1.51 

9.95e-5 

17.8 

7.0 

4.8 

3+ 

2.00 

9.95e-5 

21.4 

7.7 

5.1 

5 

0.0753 

9.95e-4 

5.2 

3.3 

2.7 

6 

0.301 

9.95e-4 

8.6 

4.7 

3.5 

7 

1.51 

9.95e-4 

17.8 

7.0 

4.8 

8 

0.0753 

3.93e-3 

5.2 

3.3 

2.7 

10 

0.301 

3.93e-3 

8.5 

4.7 

3.5 

13 

1.51 

3.93e-3 

17.4 

7.0 

4.8 

23 

0.0753 

4.03e-2 

4.3 

3.2 

2.6 

24 

0.301 

4.03e-2 

6.4 

4.3 

3.3 

25 

1.51 

4.03e-2 

9.7 

6.3 

4.6 

31 

0.301 

5.94e-2 

5.6 

4.1 

3.2 

37 

0.301 

9.86e-2 

4.3 

3.4 

2.9 

45 

0.301 

1.97e-l 

3.0 

2.7 

2.4 

47 

1.51 

1.97e-l 

4.6 

3.6 

3.2 

48 

0.753 

4.03e-l 

2.3 

2.2 

2.1 

49 

1.51 

4.03e-l 

3.0 

2.8 

2.3 

50 

3.01 

4.03e-l 

3.7 

3.0 

2.7 


Note. — We report the optimal source surface found 
by comparing pote ntial extrapola t ion an d simulation 
for all the cases of|Reville et al.| ||2015|l. The values 
vary with the stellar parameters: the rotation rate, the 
magnetic field strength and the magnetic field topol¬ 
ogy- 


field is purely radial beyond rgg, the solution of the sim¬ 
ulation has field lines that bend more gradually as the 
wind expands. The magnetic field lines can also b e colli¬ 


mated towards the rotation axis for high ro tation (Saku- 
rai|1985 Ferreira|2013 Reville et al.|2015 l, whereas this 


IS not taken into account with a potential extrapola¬ 
tion. The potential source surface model becomes inaccu- 
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It seems that globally, the potential extrapolation over¬ 
estimates the flux tube expansion, due to its inherent 
constraint to be radial beyond the source surface. At the 
pole however, the flux tube expansion is underestimated, 
at least for slow rotation. This could have consequences 
on solar wind models which derive the solar wind termi- 


rate around the optimal source surface and the deviation 

a heating with a value of 7 = 1.05 ( 

Washimi & Shibata 

from the wind solution grows with larger distances even 
though both solutions become radial. 

19931 Keppens & Goedbloed|1999l|]v 

mtt & Pudritz| 2 l)U 8 [ 

Matt et al. 

2012). The solution for 

a one dimensional. 


nal speed u sing the expansion factor (Wang fc Sheeley 


19901 

1991b 

et al. 



and references th erein) using potential ex¬ 


trapolations (see Cohen 2015|). However, as far as the 


open flux is concerned, it is always possible to find the 
optimal source surface radius that matches the simula¬ 
tion value. 

We first notice in Table that the optimal source sur¬ 
face varies with the stellar parameters (magnetic field 
strength, topology and rotation rate). Hence the fiducial 
value for the Sun, 2.hi?©, often chosen in the literature 
does not correctly predict the open flux for rapidly ro- 
tating stars with strong magnetic fields. As suggested in 
Lee et al.| (|2011|) , different values could be used f or the 
Sun, whose topology varies during one cycle (see Pinto 
et al. ( 2011 ) for a detailed study of the impact of the 11 
year solar cycle on the wind properties). 

The value of the optimal source surface radius shows 
three clear trends, which are similar to the variations 
of the average Alfven radii computed in [Reville et al.| 
(2015). First, the optimal source surface radius grows 
with the magnetic field strength. Second, it decreases 
with higher order topology. As we will see in section 
what determines the opening of the field lines is a com¬ 
petition between the magnetic forces and the thermal 
and ram pressure of the gas. The magnetic forces that 
confine the gas are proportional to the surface strength 
and follow the radial decay imposed by the topology. 
Third, we can see that rotation plays a r ole. This 
is due the magneto-centrifugal acceleration dWeber fc 


Davis|[T967| |Sakurai||1985 Mestel|p~968[ | Ustyugova et^iT 

1999[ ). S^cifically, the optimal value of rss decreases 

with higher rotation for a given vaIv^sc- The ra m pres¬ 


sure can b e significantly raised by rotation (see Reville 
et al.|2015 ) and the wind is able to open field lines closer 
to the star. 

3. PREDICTION OF THE OPEN FLUX AND 
CONSEQUENCES ON THE SOURCE SURFACE 
LOCATION 

In the previous section, we have seen that the value of 
Tas can be set such as to recover the correct amount of 
open flux. In this section we propose a method to find 
an estimate of this optimal source surface radius from 
stellar parameters. To do so we assume that a pressure 
balance is established between the flow and the magnetic 
field at the source surface and test this criteria with two 
simple wind models. 

3.1. Polytropic acceleration 


In Reville et al. (2015) we made the assumption that 
for a solar like star, the wind is driven by the pressure 
gradient of an approximately 10® K corona. We model 
this through a polytropic equation of state, mimicking 


hydrodynarnic (without magnetic field) polytropic wind 
can be computed with a Newton-Raphson method, and 
this is how we initialize our simulations. Complications 
occur when a magnetic field is introduced, especially 
with a complex topology. Semi-analytical methods have 


topologies: 

Weber & Davis 

(^1967 

) did so with one di- 

mension and later bakurai 1 


extended this result 


In our case, we want to give an estimate for the opti¬ 
mal source surface radius without having to run a MHD 
simulation. As we discussed earlier, in a wind solution 
the field lines open due to the ram and thermal pres¬ 
sure of the gas. The source surface was originally de¬ 


scribed (Schatten et al. 1969) as the radius where the 
transverse magnetic energy density becomes lower than 
the thermal energy density. From our one dimensional 
polytropic profile of the speed, density and pressure, we 
can assess the properties of the gas, and compare them 
with a “no-outflow” configuration of the magnetic field 
where none of the field lines are open (which is equiva¬ 
lent to moving rss towards infinity). We then look for an 
equilibrium between the ram and thermal pressures and 
the magnetic pressur^ 

The process is descnbed as follows. The equation: 

Phydro — F ^ram ^magt (H) 

which is equivalent to: 

P + Pv‘^ = ^. ( 10 ) 

2^0 

describes a surface in a 3D space. Our estimate is 
then simply the average spherical radius of these surface 
points. We will refer to it as the estimated source surface 
radius rss,est as opposed to the optimal source surface ra¬ 
dius rss,opt computed from our simulations. This search 
for the pressure balance is shown in Figurej^in the upper 
panel. We can see that the thermal pressure is dominant 
close to the star, but the ram pressure takes over af¬ 
ter a few stellar radii. The hydrodynamic pressure then 
crosses the magnetic pressure, hence defining the source 
surface radius estimate (although this is only a ID pro¬ 
file). As a consequence, the acceleration of the wind 
is key to determine this value, but this model does not 
yet take into account the magneto-centrifugal accelera¬ 
tion due the magnetic field lines anchored to the rotating 
star. 

3.2. Magneto-centrifugal acceleration 

The magneto-centrifugal effect is a simple cons equence 
of the exist ence of a star’s coronal magnetic field. [Schatz- 
man (1962) imagined that a magnetized wind could carry 


angular momentum, and thus could be responsible for 
main sequence sta rs braking. This concep t has been fur¬ 
ther quantified by Weber & Davis (1967) who described 


^ The magnetic pressure can be written as a tensor whose max¬ 
imum a.r nnlitiiHe onnn.serl to the p-raHien t of the magnetic field is 
iGurnett &: Bhattacharjee|2000 i. 
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Fig. 2. — Differences in the pressure comparison producing the 
estimation of the source surface radius (where the magnetic pres¬ 
sure crosses the hydro pressure) for the same v^jvesc = 1-51 and 
a dipolar topology. In the upper panel, the magneto-centrifugal 
effect is not taken into account, and the velocity profile is a so¬ 
lution of a polytropic wind. In the lower panel we used Sakurai’s 
formalism to derive the velocity profile in the equatorial plane with 
the rotation rate of case 47. We see that the acceleration is radi¬ 
cally different, hence the estimation of the source surface radius is 
much smaller when magneto-centrifugal effect is included (for fast 
rotators). 


the magnetic field as a lever-arm acting on the star’s ro¬ 
tation. Anchored to the rotating star, the magnetic field 
lines drag the gas in their rotation so that the gas feels 
a centrifugal force and is accelerated. This acceleration 
can be equivalent to, or even higher than the one due 
to the pressure gradient in thermal winds, lea ding to the 


slow and fast magnetic rotators theory (se e |Belcher & 
MacGre^^|1976^ Laiiig rs h Cassinelh] 19991. 


Weber & Uavis| (|19fc)7p used a semi-analytical method 
to compute their solution, which yields the toroidal and 
the poloidal velocity and magnetic fields near the equato¬ 
rial plane for a purely radial magnetic field. However, to 
take into account the magneto-cent rifugal effect we chose 
to implement the formalism used in Sakurai (19851. This 
formalism solves the same problem but with an improved 
methodology. We considered only the ID poloidal pro¬ 
file in this work. Details of the implementation of this 
method are given in appendix for interested readers. 



Fig. 3. — Comparison of the wind speed profiles in various cases 
of our study with a fixed magnetic field strength. The polytropic 
solution is the lowest while Sakurai Solution grows with higher 
rotation rates. Sakurai and polytropic solutions begins to differ 
with case 13 at / = 0.00393. 


In Reville et al. (2015), we showed that rotational ef- 
fects begin to be important for wind acceleration when 
/ > 0.01. To understand the strong effect of magneto¬ 
centrifugal acceleration on the velocity profile, we plot 
different wind speed solution profiles for different rota¬ 
tion rates obtained through Sakurai’s method in Figure 
As expected, the shown velocity profiles tend to the 
polytropic solution as rotation rate decreases. We can 
see that the solution begins to differ with the polytropic 
profile for case 13 (/ = 0.00393). For case 37, the ve¬ 
locity amplitude is almost four times the polytropic one 
at 10 r*. As a consequence, we expect the ram pressure 
(pram = pv^) to be strongly modified by the magneto¬ 
centrifugal effect. 

Coming back to Figure we see how the pressure bal¬ 
ance is modified when the magneto-centrifugal effect is 
included (bottom panel). We take the case of a fast mag¬ 
netic rotator (/ = 0.197, VAjvesc = 1-51) and the one 
dimensional profiles are in the equatorial plane where the 
magneto-centrifugal effect is maximum. The ram pres¬ 
sure is very significantly raised in the lower panel and 
the rss^est derived from this pressure comparison is thus 
much closer to the star. The thermal pressure is raised as 
well, due to more energy injected in the system, although 
this appears to be less significant for our pressure bal¬ 
ance. The magneto-centrifugal acceleration is also larger 
for a higher magnetic field strength, at a given rotation 
rate. 


3.3. Results 

Figure shows the estimates of and the deduced 
open fluxes for our set of parameters using both the poly¬ 
tropic and Sakurai wind solutions, and compares them 
to the outputs of the numerical simulations. For both 
panels, red points represents slow rotators (/ < 0.01) 
and for fast rotators we distinguish between the two 
wind models we implemented. Green points are com¬ 
puted with the polytropic model while blue points in¬ 
clude magneto-centrifugal acceleration through Sakurai’s 
technique. The top panel compares the optimal source 
surface radius for a given case with the source surface ra- 
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Fig. 4. — Comparison of the estimate of Vss (upper panel) and 
the corresponding open flux (lower panel) depending on whether 
the magneto-centrifugal effect is taken into account. For slow ro¬ 
tators (/ < 10“^, red points), there is no difference and we notice 
a rather good agreement between the prediction and the optimal 
value. For fast rotators, we see a large mismatch if we assume 
that the acceleration of the wind is only given by thermal gra¬ 
dient through a polytropic equation of state (green points). The 
agreement is better if we use the magneto-centrifugal wind pre¬ 
scription (blue points). The symbols stand for the topology: stars 
for dipoles, diamonds for quadrupoles, and triangles for octupoles. 
Mixed topology cases are represented with black octogones. 


dius estimate derived from formulae ([^-([l0|. We can see 
that, while both radii are close at slow rotation (red sym¬ 
bols), for faster rotation the Sakurai wind model gives 
much more accurate estimates than the polytropic one. 
For fast rotation rates and large magnetic fields, the es¬ 
timates obtained with the Sakurai wind model are three 


to four times more accurate than the estimates obtained 
by the polytropic wind model. With the Sakurai wind 
model, the average relative erroij^ of this technique is 
around 20 % (60 % with the polytropic wind model). 

The bottom panel compares the estimate of the open 
flux computed with est- It is interesting to note that 
the overall shape of this plot is inverted compared to 
the above since the open flux computed from a potential 
extrapolation decreases with increasing r^s- Once again 
the Sakurai model gives more accurate estimates. For 
large values of the open flux {^open,sim ^ 1 ) that signifi¬ 
cantly raise the magnetic torque (see formulation [^, the 
polytropic wind model almost systematically fails to re¬ 
produce the simulation value with relative error that can 
reach 160%, when we stay below 20% with the Sakurai 
model. However, pure octupolar cases with weak mag¬ 
netic fields are hard to catch with both models but they 
are unlikely to occur in realistic magnetic field configu¬ 
rations. 

Figure show some trends with the topology. For 
instance, in the upper panel, dipolar points with the 
Sakurai wind model (blue stars) are a little above the 
y = X line. In this case, the streamer is located at the 
equator where the magnetic field strength is lower and 
the magneto-centrifugal effect is maximum. Hence the 
pressure balance at the equator gives slightly better es¬ 
timates. Octupolar cases with both models (triangles) 
are usually below the y = x line. This leads to over¬ 
estimation of the open flux. It is particularly true for 
weak magnetic fields. It seems that there is a saturation 
of rss,opt around 2 i?*, i.e. the value does not go below 
this, while the pressure equilibrium can occur down to 
li?* or even yield no results if the hydro pressure is al¬ 
ways higher than the magnetic pressure. This occurs in 
our study, for all the cases with vaIv^sc = 0.0753. Those 
points are thus not represented in Figure]^ Hence to use 
this methodology, it is necessary to ensure that the mag¬ 
netic pressure is higher than the thermal pressure at the 
base of the coroneP] Also, we propose not to go below a 
saturation value of2i?* for rss,est for all cases. Doing so 
systematically improve the estimation of r^s and $open- 

Nonetheless, for complex magnetic fields, the criteria 
derived in equations ([^-p^, and illustrated in Figurej^ 
and 1 ^ works well. aHded three realistic cases that 
represents the Sun at its minimum of activity and at its 
maximum of cycle 22, and the young star TYC-5164-567- 
1. The magnetic field sp herical harmonics coeffi cients for 
the Sun are taken from |DeRosa et al. (2012) and were 
measured at the Wilcox Solar Observatory. For the Sun 
we change the value of Csjvesc to 0.26 and we consider 
a density at the surface of p* = 1.67 x 10“^® g/cm^. 
This value is calibrated such that the velocity at 1 AU 
and the mass loss rate fit observed values for 7 = 1.05, 
i.e around 450 km/s and 3 x lO“^^M 0 /yr, for both wind 
models (since the Sun is a slow rotator). We find that the 
solar Tss.opt at minimum and maximum obtained with 


ss,opt 


^ The relative error is defined by: err = 2\r, 

ss,est I / (^ ss,opt H” ^ ss,est ) 

^ This condition, equivalent to /3 < 1 at the base of the corona, 
is true for the Sun where /3 Pi 0.1. For faster r otators with higher 
magnetic fi elds, following the prescription of [Holzwarth &; Jar-| 


dine 


_ d = Pth,*/(Sj/2/Ato) OC 

(n^/Qo)^’^. Hence the condition is likely to be always true for 
the solar rotation rate and above. 
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our wind simulations bracket the fiducial value of 2.5Rq 
with rss,opt = 2.1 at maximum and rsa,opt = 3.1 at mini¬ 
mum. Our estimate at the minimum of activity perfectly 
matches the optimal value, while the estimate at the 
maximum rgs^est = is slightly under the saturation 
value of 2i?*, which has been found to be the minimum 
size of streamers in our study. 

The source surface radius is larger at solar minimum 
because of a much stronger dipole than during maximum, 
which has a strong quadrupole. Interestingly, |Lee et ah] 
(20111 predicted the opposite variation of the source sur- 
lace radius between minimum and maximum of activity. 
This latter study focused on mid-latitude coronal holes 
that are non-axisymmetric and small scale features that 
we do not account for here, to justify the variations of 
Tas- This will be investigated in the near future. 

TYC-5164-567-1 is a 120 Myr-old K-star, of mass 
M = 0.85Mq and rotational period P = 4.7 days. We 
set Csiveac = 0.285 and/9* = 4.86x10“^® g/cm^, which is 
consistent with the prescription of the coron al tempera¬ 
ture evolut i on wi th the rotation rate given by|Holzwarth 


& Jardine (2007) considering the value of 7 = 1.05 we 
use. 'i’he spherical harmonics coefficients for the surface 
magnetic field have been obtained by ZDI, using obser¬ 
vations from the spectropolarimeter ESPaDOnS (Echelle 
Spectropolarimetric Device for the Observations of Stars) 


at the CFHT (Canada Erance Hawaii Telescope) (Folsom 
et al.|[^15, submitted). Using our method and compar¬ 


ing it to a MHU simulation, we find that raa,opt = 9.8i?* 
and Taa^est = lOi?*. This example demonstrates how in¬ 
accurate the fiducial value used for the Sun can be for 
other stellar targets. This large value is mainly due to 
a 150 G axisymmetric dipole, which is common for such 
young rapidly rotating stars. The strong magnetic field 
also explains the large value of the open flux. 

For the three realistic cases, the correlation between 
the open magnetic flux and the value of Vaa is different 
than for the rest of the study. Increasing the coronal 
temperature reduces the value of r^s, at a given magnetic 
field strength due to a larger pressure gradient and more 
thermal acceleration. Those points demonstrates that 
our method is valid for different coronal temperatures. 


which are known to vary from stars to stars (^Preibisch 
1997 Gudel|[20M Holzwarth fc Jardine||2007 ). 


4. DISGUSSION 

The coronal structure of a magnetic field varies with 
stellar parameters. For a given coronal temperature, the 
magnetic streamers will grow with the intensity of the 
surface magnetic field. Magnetic topology also plays an 
important role. Hence, the PFSS model should take into 
account those parameters and we propose a method to 
do so in this paper. It might seem surprising that a one 
dimensional solution, which assumes a split-monopole 
topology, can be compared with two dimensional com¬ 
plex magnetic fields derived by 2.5D MHD wind simula¬ 
tions. This comes about because the acceleration process 
occurs on open field lines where, locally, the geometry of 
the magnetic field is close to a monopole. Hence the pro¬ 
file derived from the Sakurai technique is close to the 
one observed in simulations. This method could further 
be improved by considering a latitudinal dependency of 
the magneto-centrifugal acceleration (maximum at the 


equator) and the location of the streamers, particularly 
for 3D non-axisymmetric configurations. 

For the Sun, we find that the fiducial value of 2.5Rq 
is consistent with the optimal value we find at maximum 
and minimum of activity (2.Ii?Q and 3.17?q). However 
for younger stars with magnetic fields that can reach the 
kilogauss scale, we have seen that this value can be far 
from the optimum. The choice of has important con¬ 
sequences for the structure of the astrosphere and stellar 
dynamics. 

With physically based arguments, we propose here a 
simple way to compute the magnetic torques for any tar¬ 
get. For the Sun, using the open flux computed from 
a potential extrapolation made at the predicted by 
our technique, and using the mass loss from our wind 
solution profiles (that matches observations), we find a 
spin down time scale of 17 Gyr, at the minimum of ac¬ 
tivity. At maximum the time scale goes up to 46 Gyr. 
Those values are i n good qualitative agreement with the 
pi oneering wo r k of Skum anich J 1972) and recent s tudies 
of Matt et al. (201o|) and Gallet & Houvier (2013), from 
which a spin-down time scale ot 10 Gyr or more can be 
expected. 

For TYC-5164-567-1, we find a spin-down time scale 
of 400 Myr. Fits from observations of clusters suggest a 
value of 130 Myr ( Matt et al.|2015 ). This estimate could 
be improved by taking into account non-axisymmetric 
modes in the potential extrapolation. More complex re¬ 
constructions are also possible. In this work we only 
take into account the radial component of the magnetic 
field, but more general methods have been prop osed such 


as constant-a f orce free fields ( Her, 
potential fields (Jardine et al. 2013| 


ger [1 


1985), or non- 
his could lead 


to more accurate results tor realistic topologies obtained 
by ZDI. Moreover, our method is based on stellar param¬ 
eters that are still poorly constrained for distant stars, 
such as the density and the temperature at the base of 
the corona. The heating process used in our wind so¬ 
lution is also fairly simple and more accurate descrip¬ 
tions of the physical processes, including for instance, 
energy inputs from Alfven waves and radiative losses at 
the ba se of the corona will be implemented in the near 
futur e (Schwadro n fc McGomas|2003 Suzuki & Inutsuka 


2006 Velli 2010). 'I'he torque computation is very sen¬ 


sitive toThose prescriptions. A more detailed study of 
the torques we get with this formulation will follow in an 
upcoming paper. 

We believe this method is a step towards under¬ 
standing the coronal properties and angular momen¬ 
tum loss of low-mass stars. A open source python 
script that will perform all the calculations given a 
magnetic field strength, topology, and stellar parameter 
(i?*, M*, D*, p*, T*), can be obtained by contacting the 
first author. 
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APPENDIX 

A. WEBER AND DAVIS SOLUTION THROUGH SAKURAI TECHNIQUE 


Sakurai (19851 proposed a method to numerically compute the Weber fc Davis] ( |1967[ ) wind model, 
solves the problem of a magnetized wind anchored to a rotating star 


This method 

whose magnetic field is purely radial. The 


solution is solved in the (r, </)) plan at the equator, hence only the radial and azimuthal components of the magnetic 
and velocity fields (W, B^, i?^), and the density p and pressure p profiles are of interest. The MHD equations can 

then be integrated as follows: 


II 

■Tb 

(Al) 

= /, 

(A2) 

By = d>. 

(A3) 

(14 - nr)Br = VrB^, 

(A4) 


(A5) 

-nr)y ^ 

7 — 1 p r 2 

(A6) 


where K, /, $ and the Alfven radius, are integration constants whil e is the rotat ion rate of the star. 
Manipulating those equations can lead to a formalism, developed by Sakurai (1985), which allows one to find a 
unique solution given the stellar parameters. This appendix focuses on a description of the numerical method and the 

results, over a derivation of this formalism. _ 

Substituting all previous equations into equation |A6| we obtain an equation that only depends on r and p: 


H{r,p)=E. 


(A7) 
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Normalizing by quantities at the Alfven point, where the wind reaches the Alfven speed, we can write: 

H{r,p) = ^^H{x,y), 
rA 

where 

x = r/rA, y = p/pA, PA= 47 r/^/$^, 

P 


and 


H{x,y) = 


2x'^y‘^ 7 — 1 


0 1 w 

+ A 

X 2 


{x-l/xf _ 2 ' 

{y - 1 )^ 


P = 


AirGMpAr^ 


0 = 


w = 


lKp\ \a 

GM 

_ 


GM 


vl/ 

cy 


GM 

r 

GM 


J A 


r 

GM' 


J A 


(A8) 

(A9) 

(AlO) 

(All) 

(A12) 

(A13) 


The solution of the problem is given by the contour of H = E = E/ {GM/xa), that goes through two critical points 
{xs,ys) and {xf,yf) corresponding to the slow and fast magnetosonic points. The contour gives y{x) and thus p{r). 
The equations satisfied by H are: 


dH 

dx 

at two locations: {x,y) = {xs,ys) and {xf,yf). 



H = E, 


(A14) 


We then obtain 6 equations and 8 unknowns /3, Q,uj,E,Xs,ys,Xf,yf if we keep 7 fixed. However our unknowns are 
not independent, and are constrained by two more equations that depend on stellar parameters: 


E/i 


,1/3 = 


7 P* 
7 - Ip* 


P* \ 


Bi 


^TTPt 


GM 

r* 

7-1 


/(GMH)2/3 = 


/ 


GM\ 


2'r-4l3 


(Or*)2(4/3 


-7) 


= <?2 


(A15) 


(A16) 


The general method is thus to look for 6 parameters as a function of the two other parameters. For instance, for 
known 0 and w, a six dimensional Newton-Raphson can be used to find a unique solution. Let us define the function: 


^ fi= H{xs,ys,P) - e'^ 
/2 = H{xf,yf,P) - E 
p_ fs = da;H{xs,ys,P) 

/4 = d^H{xf,yf,P) 
fd PyB{Xg , Ps, /3) 

\ h = dyH{xf,yf,P) 

The algorithm to find the zero for this function is described as follows: 


(A17) 


• Choose an initial guess Xq 

• En = F{Xo) 

• while (Fjv > e): 

Xjv+i =Xjv- J(Xjv)-^E^ 
Ejv+i = E(Xjy+i) 
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Fig. 5.— Density profile obtained through Sakurai’s formalism in the published case of 7 = 1.2, © = 0.5, Q = 0.25 (left panel) and in the 
case 47 of our study corresponding to parameters 7 = 1.05, 0 = 1.23, 12 = 550 (right panel). In the latter case the slow and fast critical 
points are out of the domain [0,3] X [0,3]. Grey lines are the contour lines of H, the thick lines correspond to the contour at energy E 
although the black one is the only physical solution among them. 


where J is the Jacobian matrix of F taken at X^. If the initial guess is close to the solution, this method is 
remarkably fast and efficient. 

Mapping the values of E and /3 as a function of 0 and w, one can easily find the intersections of the two contour 


lines corresponding to the values qi and 92 of the resulting arrays and ^ 6 / 0 ;'^/^“'''. Then, the solution is 

fully determined from the stellar parameters, given a choice for 7 . 

In Figure]^ we show two density soluti on for different parameters, to ensure the reproducibility of our results. In 
the left panm is the solution published in Sakurai (1985), where the fast and slow critical points are close but clearly 
distinguishable from the Alfven point. The full solution for this case is: 


7 = 1.2, 0 = 0.5, a; = 0.25, /3 = 0.576, A = 1.738, (x,, ?/^) = (0.777,1.940), (x/, i//) = (1.302,0.514) 

(A18) 

In the right panel we show the density solutions for the case 47 of our study. We can see that, due to the different 
parameters, the structure of El is greatly distorted so that the critical points are out of the shown domain. The full 
solution is given by the following parameters: 


7 = 1.05, 0 = 1.235, a; = 550.0, /3 = 167.3, A = 23.0, (xs,j/s) = (0.105,902.3), (x/, i//) = (10.1,0.0066) 

(A19) 














